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ABSTRACT 

We analyze rest-frame optical morphologies and gas-phase kinematics as traced by rest-frame far- 
UV and optical spectra for a sample of 204 star forming galaxies in the redshift range 2^2 — 3 drawn 
from the Keck Baryonic Structure Survey (KBSS). We find that spectroscopic properties and gas- 
phase kinematics are closely linked to morphology: compact galaxies with semi-major axis radii r < 2 
kpc are substantially more likely than their larger counterparts to exhibit Lya in emission. Although 
Lya emission strength varies widely within galaxies of a given morphological type, all twelve galaxies 
with Lya equivalent width WLya > 30 A have small semi-major axes. The velocity structure of 
absorption lines in the galactic continuum spectra also varies as a function of morphology. Galaxies 
of all morphological types drive similarly strong outflows (as traced by the blue wing of interstellar 
absorption line features), but the outflows of larger galaxies are less highly ionized and exhibit larger 
optical depth at the systemic redshift that may correspond to a decreasing efficiency of feedback in 
evacuating gas from the galaxy. This u ~ km s^^ gas is responsible both for shifting the mean 
absorption line redshift and attenuating W^Lya (via a longer resonant scattering path) in galaxies with 
larger rest-optical half light radii. In contrast to galaxies at lower redshifts, there is no evidence for 
a correlation between outfiow velocity and inclination, suggesting that outflows from these puffy and 
irregular systems may be poorly collimated. Our observations arc broadly consistent with theoretical 
models of inside-out growth of galaxies in the young universe, in which typical z ^ 2 — 3 star forming 
galaxies are predominantly unstable, dispersion-dominated, systems fueled by rapid gas accretion 
that later form extended rotationally-supported disks when stabilized by a sufficiently massive stellar 
component. 

Subject headings: galaxies: high-redshift — galaxies: fundamental parameters — galaxies: structure 



1. INTRODUCTION 

Galaxy evolution is a process driven largely by the 
lifecycle of gas and its accretion, conversion into stars, 
and eventual loss via supernova-driven winds. Dur- 
ing the peak epoch of galaxy formation {z ^ 2 — 3 ; 
IDickinson et al.lboOStlReddv et al.l20"ollZhu et al.ll2009D 
those systems that dominate the luminosity function are 
especially gas-rich, with gas frequently acco unting for 
50% o r more of the total b aryonic mass fe.g.. lErb et al.l 
I2006cl: iTacconi et al|[2010l ). This large supply of hy- 
drogen fuels sustained sta r formation rates (SFR) com- 
monl y ~ 30 Mq yr-i fe.g.. lErb et al.ll2006bHWuvts et al.l 
1201 ID and SFR surface densities pe aking above ~ 1 Mq 
yr~^ kpc~^ fe.g.. lGenzel et aLll2011[ ) that are comparable 
to those in transi ent starbursts in the modern universe 
(|Kennicuttill998bD . 
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The mechanisms by which these galaxies acquire such 
large quantities of gas are uncertain; although gas has 
traditionally b een assumed to accrete primarily through 
mergers (e.g., [Robertson et al.l 120061) and 'hot mode ' 
accretion via yirial shocks Ce.g.. iRees fc Ostrikeii 119771 : 
IWhite fc Reei [19781: IWhite fc FrenkI I1991D . some re- 
cent arguments suggest that semi-cont inuous cold accre- 



tion f r om cosmological fi l aments (e.g.JKeres et al.l 
20091: IDekel et al.l [20091 [Bournaud fc Elmegreed 



2005 



2009; 



Ceverino et al. n2oioir may be required to sustain the 



observed SFR and produce the structures observed at 
z ~ 2. As discussed by Reddy et al. (2012) however, 
star formation at 2; ^ 2 — 3 appears to have been an in- 
efficient process, indicating that SFR is unlikely to have 
been limited primarily by the cold gas accretion rate. 

In contrast to galaxies in the nearby universe, z ^ 2 
star forming galaxies do not exhibit r egular spiral struc - 
ture (with some notable exceptions; [Law et all [2012bO 
and instead tend to be compact, triaxial systems con- 
sisting of o ne or more irregularly - shaped clumps of emis- 
sion (e.g.. [Conselice et al.l l2005[: [Elmegreen et al.[ [2005[ : 
[Ravindranath et al.l [2006[ : [Law et al.l [2012al : and refer- 
ences therein). In many of the most massive galaxies 
(M, > 3 X IO^^Mq) such clumps may be embedded in low 
surface-brightness structures reminiscent of thick disks 
with typical scaleheights ^ 1 kpc (e.g., Elmegreen 
fc Elmegreen 2006; Genzel et al. 2006, 2008; Forster 
Schreiber et al. 2011a, 2011b) sustained by their large 
vertical velocity dispersions (e.g., Genzel et al. 2006; 
Law et al. 2007, 2009; Forster Schreiber et al. 2009; 
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Swinbank et al. 2011). 

The large cold gas fractions of these galaxies appears 
to play a significant role in their dynamical and mor- 
phological evolution. Simulations (e.g., Noguchi 1999; 
Immeli et al. 2004; Bournaud et al. 2007; Bournaud 
& Elmegreen 2009) indicate that such gas-rich galax- 
ies should be highly unstable to fragmentation, form- 
ing clumps similar to those observed which may either 
migrate to form a central bulge (Immeli et al. 2004; 
Elmegreen et al. 2008; Ceverino et al. 2010) or be rapidly 
disrupted by stellar feedback (Wuyts et al. 2012). A 
reliable understanding of the relation between galactic 
stellar morphology and its gas-phase properties is there- 
fore critical for constraining models of galaxy formation 
in the young universe. 

While observational signatures of infalling gas 
(whether hot or cold-mode) have proven subtle and dif- 
ficult to detect (although see lKulas et al.l[2012[ ). the in- 
tense star formation of the galaxies gives rise to ubiqui- 
tous and well-studied galactic-scale outflows evident in 
rest-frame UV spectra (e.g., Shapley et al. 2003; Weiner 
et al. 2009; Steidel et al. 2010). These outflows are 
frequently powerful enough to drive enriched gas to the 
virial radius and beyond with mass loss r ates compara- 
ble to or exceed i ng the star formati on rate (jPettini et al.l 
[2000t l&bl [20081: iSteidel et al.ll2010D . Numerous authors 
(e.g., Shapley et al. 2003, Penterrici et al. 2010; Ko- 
rnei et al. 2010; Steidel et al. 2010; Jones et al. 2012) 
have therefore sought to use the information about op- 
tical depth, covering fraction, geometry, and gas kine- 
matics encoded within such spectra to learn about the 
mechanisms by which cold gas is coupled to feedback and 
physical observables such as stellar mass, star formation 
rate, and stellar population age. Efforts to understand 
the role of the gas in shaping the evolution of galaxy 
morphology however (e.g.. Law et al. 2007b; Penterrici 
et al. 2010) have historically been limited in their ability 
to test current galaxy formation models by the lack of 
availability of high resolution rest-frame optical imaging 
tracing the bulk of the stellar mass. 

In the present contribution, we combine rest-frame op- 
tical imaging of a large sample of z ^ 2 star forming 
galaxies recently obtained with the WFC3 camera on 
board the Hubble Space Telescope (HST) with ground- 
based rest-UV spectroscopy to investigate the relation 
between galaxy morphology and gas-phase kinematics. 
In Section [2] we describe the target galaxy sample and 
detail the methods by which morphological, photomet- 
ric, and spectroscopic data were obtained. We discuss 
the Lya emission properties of the sample in and the 
structure of major interstellar absorption lines arising in 
outflowing cold gas in fJH highlighting the changing ve- 
locity substructure of the absorption lines as a function 
of galaxy morphology. We discuss the implications of our 
findings for the evolution of typical z ^ 2 — 3 star forming 
galaxies in S}5l combining our results with recent obser- 
vations from IFU kinematic surveys. We summarize our 
conclusions in [J6l 

Throughout our analysis, we adopt a standard ACDM 
cosmology based on the seven-year WMAP results (Ko- 
matsu et al. 2011) in which Hq = 70.4 km s~^ Mpc"\ 
= 0.272, and = 0.728. 

2. OBSERVATIONAL DATA 



Our galaxy sample is drawn from the Keck Baryonic 
Structure Survey (KBSS; Trainor & Steidel 2012), a cat- 
alog of z ~ 1.5 — 3.6 star forming galaxies selected 
according to optical UnGTZ color fStcidcl et al. 200^ 
120041 lAdelber ger et al.ll2004D and confirmed using rest- 
UV spectroscopy. These galaxies typically lie near sight- 
lines to hyperluminous background QSOs (zqgo ~ 2.7) 
scattered widely across the sky. We define two specific 
samples of galaxies used in our analysis: 

Parent sample:: A sample of 204 galaxies with 
HST/^YCi imaging data, magnitudes iJieo < 24 
AB, and redshifts z ~ 1.5 — 3.5 (see Fig. [T]) derived 
from rest-UV spectroscopy. 

HAHQ subsample:: A subgroup of 35 galaxies from 
the parent sample for which the rest-UV spec- 
troscopy is particularly high-quality ('uvqual' > 1; 
see §2.2|) and secure systemic redshifts in the range 
z = 2.0 — 2.5 (see Fig. [T]) have been obtained from 
Ha nebular emission line spectroscopy. 

2.1. Rest-Optical Morphologies 

As part of Hubble Space Telescope {HST) Cycle 17 pro- 
gram GO-11694, we obtained WFC3/IR imaging of 306 
UnGTZ color-selected star forming galaxies with spectro- 
scopic redshifts in the range z = 1.5 — 3.6 drawn from the 
KBSS. The details of these observations have been de- 
scribed at length by Law et al. (2012a; hereafter 'Paper 
I'). In brief, we used the F160W (Acff = 15369 A) filter 
to trace rest-frame optical emission (A ^ 3400 — 6100 A 
for the parent sample, A - 4400 - 5100 A for the HAHQ 
sample). Our data reach a depth of 27.9 AB for a 5a de- 
tection within a 0.2 arcsec radius aperture (representing 
a total of 8100s integration in each pointing), and have 
a PSF FWHM of 0.18" that corresponds to 1.5 kpc at 
z ~ 2. 

The morphologies consist of three qualitative types. As 
defined in Paper I, these are 

Type I:: Single nucleated component with little to no 
low surface brightness emission. 

Type 2:: Two or more distinct nucleated components of 
comparable magnitude and little to no low surface 
brightness emission. 

Type 3:: Extended objects with non-axisymmetric low 
surface-brightness features. 

Numerous morphological statistics have been defined 
in the literature that attempt to encapsulate this range 
of morphologies quantitatively. In Paper I we calculated 
7 of the most widely used such statistics for each of the 
galaxies in our WFC3 imaging fields. These statistics 
are: 

r:: Semi-major axis radius derived from fitting Sersic 
(1963) models convolved with the observational 
PSF to the F160W hght profile using the GAL- 
FIT (Peng et al. 2002, 2010) analysis package. 
We also defi ne the circularized effective radius 
(re = r^jbja) for consistency with previous anal- 
yses, where bja is the apparent minor/major axis 
ratio. 
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Fig. 1. — Redshifts, stellar population parameters, and effective radii derived for the parent (black open boxes) and Ho high-quality 
(HAHQ; red filled boxes) samples of 204 and 35 galaxies respectively. The two samples have significant overlap, but the HAHQ sample is 
more restricted in redshift range, and has marginally larger stellar masses and population ages than the parent population on average. 



n:: Radial index of the best-fit Sersic model profile. 

G:: Gini statistic (Abraham et al. 2003; Lotz et al. 2004) 
quantifying the cumulative flux distribution among 
the galaxy pixels. 

Af2o:: Second order moment of the spatial distribution 
of the light profile (Lotz et al. 2004, 2006). 

C:: Concentration index of the flux distribution about 
the center of the galaxy (Bershady et al. 2000; 
Conselice 2003). 

A:: Asymmetry index (Schade et al. f995; Conselice et 
al. 2000) quantifying the 180° rotational asymme- 
try of the galaxy. 

5':: Multiplicity index (Law et al. 2007b) quantifying 
the 'work' required to assemble the distribution of 



pixel fluxes within a galaxy (distinguishing single- 
component from irregular and multiple-component 
sources). 

Although the first two morphological statistics (r, n) 
assume a parametric model convolved with the obser- 
vational PSF, the latter five (G, M20, C, ^, ^') are non- 
parametrically derived from the spatial arrangement of 
individual pixel fluxes. The values calculated for these 
five statistics depend systematically on the adopted seg- 
mentation map defining which pixels are considered to 
be part of a given galaxy. Our approach is placed on 
the systematic reference frame of Lotz et al. (2008a, 
2008b, and references therein) using the transformations 
described in Appendix A.l of Paper L 

All seven morphological statistics are robust for galax- 
ies with HiQo < 24 AB, and show systematic bias at 
fainter magnitudes (see Appendix A. 2 of Paper I). We 
therefore restrict our attention in the present contribu- 
tion to those 204 galaxies with i/ieo < 24.0 and no ev- 
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idence of QSO or AGN in their UV spectra or 
broadband photometry fi )2.4p . 

2.2. Rest-UV Spectroscopy 

Rest-frame UV spectra for each of our galaxies have 
been obtained using the LRIS multisht spectrograph on 
Keck I. Details concerning the observational setup, data 
reduction, and calibration process of these spectra have 
been presented at length elsewhere in the literature (e.g., 
Shapley et al. 2003; Steidel et al. 2010; and refer- 
ences therein). In general, individual spectra represent 
integration times ^1.5 hours with characteristic veloc- 
ity resolution ares ~ 160 km s~^ . These spectra vary 
considerably in quality according to both the intrinsic 
brightness of the target galaxy and the conditions under 
which particular slit masks were observed. We sort the 
galaxies according to the quality of their spectra, assign- 
ing them quality indices 'uvqual' ranging from (one or 
two spectroscopic features identified with relatively low 
confidence) to 3 (multiple strong spectroscopic features 
identified with extremely high confidence). 

The limited S/N ratio of even the highest-quality spec- 
tra however means that it is typically only possible to 
estimate the mean interstellar absorption (emission) line 
redshift zjg (zem) for individual galaxies. In order to in- 
vestigate the detailed structure of spectral lines it is nec- 
essary to stack together multiple spectra to reach higher 
S/N ratios. This stacking process is described in detail 
by Steidel et al. (2010), and includes doppler-correcting 
each spectrum to the rest framed resampling to a com- 
mon 0.33 A/pixel spectral dispersion, and rescaling to a 
common mode in the wavelength range A ^ 1250 — 1500 
A before combining individual spectra using a sigma- 
clipped mean algorithm. 

In Figure [5] we show the composite spectrum for the 
204 galaxies in our parent sample. This composite spec- 
trum exhibits a mixture of Lya absorption and emission, 
deep blended interstellar absorption features due to var- 
ious metal species in the galactic ISM, and a host of 
weaker features due to stellar photospheric absorption 
and fine structure emission. These absorption features 
are typically blueshifted with respect to the systemic red- 
shift (corresponding to absorption arising in the near side 
of an outfiowing galactic wind), while Lya emission is 
redshifted (corresponding to off-resonance backscattering 
from the far side of the outflow) . As described at length 
by Steidel et al. (2010; see also Kulas et al. 2012), this 
outflow is likely composed of discrete gas clouds whose 
expansion velocities increase with galactocentric radius. 

2.3. Rest-optical Spectroscopy 

Although the rest-UV spectra give an indication of the 
redshift of the interstellar absorption lines zjg, it is not 
possible to perform a detailed analysis of the kinematic 
structure of the outflowing gas without knowing the sys- 
temic redshift Zgya of the galaxy as well. We therefore 
define the 'Ha - high quality' (HAHQ) subsample of 
galaxies with particularly robust measurements of zis 
and for which accurate (uncertainty ^ 60 km s~^ ; Stei- 
del et al. 2010) systemic redshifts have been derived 

Where available we use the nebular emission li ne r edshift de- 
rived from rest-optical spectroscopy as detailed in il2.3l otherwise 
we adopt the prescriptions given by Steidel et al. (2010). 



from Ha nebular emission using a combination of long- 
slit Keck/NIRSPEC (Erb et al. 2006b), integral-field 
Keck/OSIRIS (Law et al. 2007a, 2009; Wright et al. 
2009) and VLT/SINFONI (Forster Schreiber et al. 2009) 
spectroscopy. The HAHQ sample represents the subset 
of 35 of the 89 galaxies considered by Steidel et al. (2010) 
for which _ffS'T/WFC3 morphologies have been obtained. 

2.4. Stellar Population Models 

Extensive KBSS ancillary data exists for our target 
galaxies, typically including deep ground-based optical 
and infrared J/ Kg photometry, and in some cases Spitzer 
IRAC and/or MIPS photometry. In combination with 
HiQo magnitudes derived from the HST/WFC3 imaging 
data, the long- wavelength photometry permits us to con- 
struct stellar population models by fitting the broadband 
spectral energy distribution of each galaxy as detailed 
in Paper I (see also Shapley et al. 2001, 2005; Erb et 
al. 2006c; Reddy et al. 2006, 2010). In brief, we fit 
Chariot & Bruzual (2012, in prep.) stellar population 
synthesis models in combination with a Chabrier (2003) 
IMF, Calzetti et al. (2000) extinction law, and a con- 
stant (r = oo) star formation history using a customized 
IDL code (Reddy et al. 2012). 

As illustrated in Figure [1] galaxies in the parent 
sample have stellar masses in the range ^ 10^ — 
10^^ M0, star formation rates ^ 30AfQ yr~^, and 
moderate color excess E{B — V) ^ 0.0 — 0.4. The 
HAHQ subsample is biased relative to the parent sam- 
ple by the requirement of both detecting Ha emission 
and observing a high-quality UV spectrum (favoring 
UV- bright objects), and therefore has a slightly larger 
mean stellar mass ((log Af*/MQ)HAHQ = 10.1 ± 0.4 vs. 
(log M./Mq) 

Parent — 

9.9 ±0.50 

3. LYa EMISSION 

Lya emission originates deep within individual star 
forming regions and is resonantly scattered throughout 
the interstellar medium until either being absorbed on 
a dust grain or redshifted sufficiently far off-resonance 
that it can escape the galaxy. Although almost all star 
forming galaxies app ear to have Lya em itting haloes on 
large angular scales (jSteidel et al.ll201ll ). the Lya profile 
of the central few kpc traced by slit spectra is dependent 
upon the structure, kinematics, and ionization proper- 
ties of each galaxy and its interstellar and circumgalac- 
tic medium. Substantial efforts have been made over the 
last decade to attempt to disentangle these various effects 
and connect the central Ly a line profile w i th glo bal stel- 
lar population properties. IShapley et al.l (|2003l ). for in- 
stance, found that z ~ 3 Lyman Break Galaxies (LBGs) 
with stronger Lya emission tend to have lower SFR, 
bluer UV continua, and weaker low-ionization interstel- 
lar absorption lines. Similar results were foun d for z ^ 2 
star for ming galaxies by Erb et al.l (|2006a[ ) and lLaw et al.l 
(|2007b( l. who noted that galaxies with stellar masses 
M^, > 1O^°M0 on average have weaker Lya emission 
than do galaxies with M^, < lO^^M©, and at z ~ 4 by 
Jones et al. (2012). These results were borne out by 
the more detailed analysis of iSteidel et al.l (|2010( ). which 
demonstrated that Lya emission is both weaker and more 

^ la values represent standard deviation about the mean. 
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Fig. 2. — Flux-calibrated rest-UV Keck/LRIS-B spectrum constructed from stacking 204 star-forming galaxies in our parent sample in the 
redshift range 1.5 < z < 3.5. Vertical dotted lines indicate the fiducial wavelengths of major emission and/or absorption line transitions. 



strongly redshifted from the systemic velocity in galax- 
ies with high baryonic mass, consistent with the picture 
that Lya emission from the deepest gravitational poten- 
tial wells must scatter further from resonance before be- 
ing able to escape. Using a similarly selec t ed sa mple of 
LBGs at z = 2.5 - 3.5, IPentericci et~all pOTot ) found 
comparable results, showing that galaxies without Lya 
emission tend to be more massive and dustier, although 
these authors observed that stellar population ages and 
SFRs did not depend strongly on Lya emission charac- 
teristics. 

The relation between Lya emission and rest-frame 
UV gal axy m orphology has previousl y been i nvest i- 
gated bvlLaw e t al. (20073) and Pentcri cci et all (|2010D . 
iLaw et al.l ()200 7b) found that Lya emission was corre- 
lated with the Gini coefficient in the sense that high-G 
(i.e., strongly nucleated) systems were more likely to ex- 
hibit Lya emission, which we previously ascribed to the 
lower q uantity of dust in strong ly nucleated galaxies. Al- 
though [Pe^ericcr^F^L| (|201 0') found no significant cor- 
relation of Lya emission with any of the non-parametric 
morphological indices {C,G,A, or M20), these authors 
noted that while galaxies with weak or no Lya emission 
spanned a range of projected angular sizes, galaxies with 
strong Lya emission tended to be exclusively small. 

In ij3.1l we discuss the probability of galaxies within our 
sample to exhibit Lya emission as a function of various 
physical properties, and expand our discussion in tJ3.2l to 
consider the trends in emission line strength within the 
Lya-emitting subset. 

3.L Probability of Lya Emission 

We divide the parent sample into two groups accord- 
ing to whether Lya emission is present (71 galaxies) or 
not (117 galaxies) in the UV spectrum of each galaxyQ 
In Table [T] we list probabilities calculated using the non- 

* Note that this is not a selection in equivalent width akin to that 
used by, e.g., Gawiser et al. (2007). Rather, our sample include 
galaxies with Lyo emission of any strength, including those with 
a small emission peak superimposed on a large Lya absorption 
trough. We rejected 16 galaxies in the redshift range 2 = 1.5 — 1.7 
from our sample since Lya emission at these redshifts would not 
fall within the LRIS UV spectral range. 



parametric Kolmogorov-Smirnov (KS) test to evaluate 
whether the two groups are statistically consistent with 
the null hypotheses that they are drawn from the same 
distribution of stellar mass, SFR, morphology, etc. We 
find that SFR and the morphological parameters M20, 
C, A, Vt, and n are poor predictors of the likelihood 
of observing Lya emission from a galaxy; galaxies with 
and without Lya in emission have statistically indistin- 
guishable distributions. In contrast, the null hypothe- 
sis is rejected at greater than 94-95% confidence when 
we examine the galaxy mass, age, color excess, G, or r 
distributions for Lya emitters and non-emitters. These 
correlations are in the sense that galaxies with Lya emis- 
sion on average are smaller, have slightly lower mass 
((log(Af*/AfQ)) = 9.82 vs 9.96), are slightly younger 
(500 vs 800 Myr), and less dusty {{E{B - V)) = 0.18 vs 
0.23) than galaxies with no Lya emission. These results 
are broadly consisten t with relation s that h a ve already 
been est ablished b y rShapley et all (I2003D lErb et al.l 
(2006 31). lLawet al.1 (l2007bn . ISteidel erarr(|2010D . and 

Pentericci et al.l (|2010DT 

We note, however, that iKornei et al.l ()2010l ) studied 
the Lya properties of z 3 LBGs and found a conflict- 
ing result that z ~ 3 Lya emitters tended to have stellar 
populations slightly older than their non-emitting kin, 
which was hypothesized to indicate that Lya emitters 
may represent a later stage in galaxy evolution in which 
supernova-driven outflows have reduced the dust cover- 
ing fraction (see also Shapley et al. 2001). The difi^erence 
with our present findings may partly reflect the different 
mean redshift range of the samples and the manner in 
which Lya emitters are defined; although we define a 
Lya emitter as a galaxy with any observab l e emis sion 
component to the Lya profile, IKornei et al.l (|2010( ) in- 
stead used an equivalent width measure that is sensi- 
tive to both the relative strength of Lya emission and 
absorption components. Alternatively, the relation be- 
tween Lya emission strength and stellar population age 
may simply be too weak (96% rejectio n of the null hy- 
pothe sis) to measure consistently; while [Pentericci et alJ 
(|2010( ) concurred that Lya emitters tended to be slightly 
younger than similar-selected non-emitters on average, 
these authors too noted that the correlation was only 
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TABLE 1 

KS Probability that Lycv emitters and 
non-emitters are drawn from the same 
distribution^ 



Quantity Pks 



M. 6% 

Age 8% 

SFR d 99 % 

SigmaSFR " 1 % 

E(B-V) * 3% 

G 8% 

M20 83% 

C 26% 

A 99% 

* 56 % 

r g 3% 

re 0.5% 

n 15 % 

6/a ' 59 % 

Hieo 12 % 



^ All values are the probability of the null hypothesis that 
the quantities are drawn from a statistically identical dis- 
tribution of the indicated parameter. The sample size is 
188 galaxies. 

^ Stellar mass, estimated from SED fitting. 

^ Average population age from a constant SFR SED 

model. 

^ Star formation rate measured from SED fitting. 

^ Average star formation rate surface density Ssfr — 

^ Color excess due to dust extinction from a constant SFR 
SED model. 

g GALFIT semi-major axis radius. 

^ GALFIT circularized effective radius. 

^ Morphological minor/major axis ratio. 

marginally significant. 

As indicated by Table [1] within our sample of star 
forming galaxies there is no statistically significant cor- 
relation between the presence/absence of Lya emission 
and the morphological indi ces C, A, M^oi ^1 o r n. Simi- 
lar to the result obtained bv lLaw et al.l ()2007bl ) , there is a 
relation (92% confidence) between Lya emission and the 
Gini parameter G in the sense that Lya emitters tend to 
have marginally more strongly nucleated light profiles in 
the rest-frame optical. This trend is far from compelling 
however and the mean values of G for galaxies with and 
without Lya emission respectively are only minimally 
different ( (G)Lya = 0.49 vs (G')NoLya = 0.48). 

In contrast, the presence of Lya emission is strongly 
correlated with the morphological semi-major axis ra- 
dius r, circularized radius re, and (by construction) the 
star formation rate surface density Esfr = SFR/7rr^. 
These correlations are significant at the 97-99% confi- 
dence level, comparable to or stronger than the well- 
known anti-correlation between Lya emission and dust 
content E(B — V) ). Similarly to th e rest- U V angular size 
relation noted by iPentericci et al.l ()2010[ ) , Lya emitters 
are more likely than non-emitters to have small sizes, 
and the smallest galaxies in our sample are much more 
likely to have Lya in emission than are the largest galax- 
ies with r > 3 kpc (see Figure [3|). This is unlikely to be 
due to any intrinsic correlation between semi-major axis 
radius and dust content; although r and E{B — V) of 
the parent sample are marginally correlated at the '~ 2a 
level (based on the Spearman rank correlation test), this 



correlation is in the sense that larger galaxies tend to be 
marginally less dusty. 

We illustrate the spectroscopic difference between 
large- and small-r galaxies by stacking the halves of the 
parent and HAHQ samples with the smallest and largest 
semi-major axis radii (Figure |4])I3 The Lya emission 
equivalent width H^Lya of the large and small-r stacks is 
1.3 A and 2.6 A respectively in the parent sample, repre- 
senting a factor ^ 2 increase in the mean Lya emission 
strength^ This differ ence is even more pronounced in 
the HAHQ subsamplc, increasing by a factor ^ 10 from 
the large-r (M^Lya = 0.7A) to small-r (Wiya = 6.7 A) 
stacks, corresponding to a K-S probability of ~ 0.2% that 
the radii of Lya emitters and non-emitters are drawn 
from the same distribution. 

3.2. Strength of Lya Emission 

Within the Lya-emitting subsample similar correla- 
tions exist between Lya emission strength and physical 
galaxy properties. We compute the emission line equiv- 
alent width Wi^ya relative to the continuum level in our 
flux-calibrated spectra for each of the 71 galaxies with 
such emission in the parent sample. In Table [2] we tabu- 
late the correlations between M^Lya and physical galaxy 
properties in terms of the number of standard devia- 
tions from the null hypothesis (that the two variables 
are uncorrelated) using the Spearman rank correlation 
test. The factors most strongly influencing W^Lya are the 
dust content E{B — V) of a, galaxy (2.4cr), its total star 
formation rate (2.4cr)I3 and its physical size (2.3cr), in 
the sense that larger, dustier galaxies with higher SFR 
have lower Whya- 

We illustrate these three trends in Figure [H Similarly 
to previous studies (e.g., Shapley et al. 2003; Kornei et 
al. 2010; Penterrici et al. 2010), we find that galaxies 
with the strongest WLya tend to contain lower than av- 
erage quantities of dust and have lower star formation 
rates. Additionally, galaxies with large rest-optical sizes 
(r > 3 kpc) tend to have W^ya < 30 A, while the smallest 
galaxies (r < 2 kpc) exhibit a much wider range in equiv- 
alent width from Whya = — 80 A. Star-forming galaxies 
in our optical-color selected sample with Wi,ya > 30 A all 
have small rest-optical radii and exclusively Type I and 
II morphologies (i.e., one or two nucleated components 
to the light profile with little to no extended and/or dif- 
fuse low surface brightness emission in the WFC3 imag- 
ing data; see Figure [6]). The correlation between Lya 
equivalent width and optical semi-major axis radius is 
independent of the correlation with E{B — V); r and 
E{B — V) are consistent with the null hypothesis of be- 
ing uncorrelated (O.lcr rejection of the null hypothesis 
using the Spearman rank correlation test), and confin- 
ing our sample to a narrow range of values in E{B ^ V) 

^ The median radius used to divide the samples is 2.08 kpc and 
1.98 kpc for the parent and HAHQ samples respectively. 

By comparison, Monte Carlo tests stacking randomly selected 
halves of the galaxy sample indicate that H^Lya of the stacks typ- 
ically vary by ~ 30% for the parent sample, and by a factor ~ 2 
for the HAHQ subsample. 

Although the SFR has no statistically significant relation to 
the probability of observing Lya in emission, for the subsample of 
galaxies with Lya in emission we find the SFR to be correlated 
with the strength WLya of the emission. 
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TABLE 2 

Correlation of Wi^ya and Galaxy Properties^ 



Quantity 


Std. Dev. 




-1.7 


Age " 


0.5 


SFR <^ 


-2.4 


SSFR 


-0.2 


E(B-V) f 


-2.4 


G 


-0.2 


M20 


2.2 


C 


-0.7 


A 


1.4 




1.6 




-2.1 


r h 
' e 


-2.3 


n 


-1.1 


b/a ' 


-1.9 




1.4 



Fig. 3. — Left panel: Histogram of semi-major axis radius r for galaxies with and without Lya emission component in the 183 galaxies 
of the parent sample for which the LRIS spectral range includes the Lya transition. Fractions shown are relative to the total number of 
galaxies with (71 galaxies) and without (117 galaxies) Lyo emission respectively. Right panel: Fraction of galaxies in each radial bin that 
show Lya emission; Icr uncertainties are estimated using a binomial probability distribution. 

sion, galaxies with a measurable Lya emission compo- 
nent tend to be physically smaller systems with lower 
stellar mass and less dust, with marginally younger stel- 
lar populations (although c.f. Kornei et al. 2010 for 
galaxies at z ^ 3), lower star formation rates, and a 
slightly greater degree of nucleation in the rest-frame UV 
and optical light distribution. 

3.3. Galaxies with Exclusively Lya Emission 

Most commonly, galaxies with visible Lya emission 
also show some absorption features in their rest-UV spec- 
tra as well (albeit frequently weak; see discussion by 
Shapley et al. 2003). Of the 71 galaxies in the parent 
sample with Lya emission, 7 do not have any identifiable 
absorption features (even in a composite stack). While 
the Lya emission equivalent width for these sources range 
over a factor - 20 (WLya = 17 - 184 A), ah seven have 
Type II morphologies composed of two or more widely 
separated clumps (see Figure |6]) , accounting for nearly 
half of the 17 total galaxies in the parent sample with 
Lya emission and Type II morphology. By comparison, 
100% of the 28 (26) galaxies with Type I (III) morphol- 
ogy and Lya emission have visible absorption features in 
their UV spectra, indicating that multi-component mor- 
phology is a strong factor in determining whether absorp- 
tion line features will be observed in rest-UV spectraQ 
Morphological indices measuring the two-dimensional 
distribution of the light from a galaxy (e.g.. A, M20, 
^f) therefore easily identify galaxies with exclusively Lya 
emission in their rest-UV spectra and suggest that these 
multi-component galaxies are prime major-merger candi- 
dates. If these galaxies are indeed mergers caught in the 
nearby-pair phase, the disturbed ISM of the combined 
system may result in a reduced gas covering fraction, 
making it easier for resonantly-scattered Lya emission 

The fraction of galaxies with Type II morphology and no Lyo 
emission that do not have visible rest-UV absorption features is 
not possible to determine, since by definition such galaxies would 
not have spectroscopically confirmed redshifts. 



^ All values arc the number of standard deviations from 
the null hypothesis that the quantities are uncorrelated, 
based on a Spearman rank correlation test. Positive (neg- 
ative) values indicate positive (negative) correlation be- 
tween the quantities (i.e., larger r corresponds to lower 
^^Lyoi)- The sample size is 71 galaxies. 
^ Stellar mass, estimated from SED fitting. 
^ Average population age from a constant SFR SED 
model. 

^ Star formation rate measured from SED fitting. 

^ Average star formation rate surface density Ssfr — 

SFI1./nr^. 

^ Color excess due to dust extinction from a constant SFR 
SED model. 

S GALFIT semi-major axis radius. 

^ GALFIT circularized effective radius. 



^ Morphological minor /major axis ratio. 

results in a similarly strong trend of Whya with r after 
accounting for the reduced sample size. 

Combining our results with those previously published 
in the literature, the broad picture of z 2 — 3 Lya 
emitters is therefore relatively clear. Compared to star 
forming galaxies at this epoch with negligible Lya emis- 
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Fig. 4. — Continuum-normalized stacked Keck/LRIS-B UV spectra of the 50% of galaxies with smallest (blue lines) and largest (red 
lines) sizes for the parent (200 galaxies) and HAHQ (35 galaxies) galaxy samples (note differing vertical scales). The vertical dotted lines 
indicate the wavelengths of major spectral features arising from a combination of low/high ionization interstellar absorption and various 
emission mechanisms. While not labeled, the emission feature at ~ 1300 A is the fine-structure emission line Si II*A1309A. Note the drastic 
difference in Lya profiles and varying depths of the low-ionization absorption features. 



to escape from the system. However, type II galaxies 
are no intrinsically brighter in Lya emission than Type 
I or III galaxies; instead, the only significant difference 
is that type II galaxies are less likely to have absorption 
lines in their UV spectra. A simpler explanation may be 
that when galaxies with multiple components separated 
by distances ^ 1 arcsec are observed with ground-based 
spectrographs the slit tends to be placed midway between 
the two components, thereby systematically reducing the 
S/N ratio of the continuum emission against which ab- 
sorption features are measured. As illustrated in Figure 
[6l in 50% of cases the majority of the rest-UV contin- 
uum flux is expected to fall outside the slit based on the 
position angles used for the LRIS spectroscopy; even in 
cases for which the slit was aligned approximately along 
the separation vector (e.g., Q1623-BX547) the greater ef- 
fective area over which the light is distributed may still 



increase the noise in the continuum such that the obser- 
vations are systematically biased against the detection of 
absorption lines. 

4. METAL LINE FEATURES 

The deepest absorption lines in the rest-UV spectra 
are the Sin, Oi, Cii, Fell, and Alii lines that arise in 
the metal-enriched neutral ISM, which typically trace 
massive galactic-scale outflows driven by the intensely 
starbursting galaxy. As noted by previous authors (e.g., 
iShaplev et al. 2003) there is generally an anti-correlation 
between Lya emission and interstellar absorption line 
strength; since the metal transitions are optically thick 
declining absorption line strength indicates a declining 
gas velocity dispersion and/or covering fraction, which 
also permits a greater Lya escape fraction. 

In order to understand the physical mechanisms un- 
derlying variations in absorption line strength in greater 
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Fig. 5. — Lya emission equivalent width VKLya as a function of rest-optical semi-major axis radius r, dust content E[B — V), and 
SED-derived star formation rate for the 71 galaxies with measurable Lya emission in the parent sample (black points). Red points denote 
the 12 Lyo emitting galaxies that are also in the HAHQ subsample. Wi^ya decreases with r, E{B — V), and SFR; the significance of the 
correlations as estimated from the Spearman rank correlation test are indicated in each panel. 
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Fig. 6. — HST/WFC3 F160W rest-optical morphologies of the twelve galaxies in the parent sample with Lyo emission equivalent width 
Wi^ya > 30 A (top panel) , and of the 7 galaxies in the parent sample with Lya emission but no absorption features in their rest-U V spectra. 
Both sets of galaxies are sorted in order of decreasing WLya)- Note that all twelve galaxies in the top panel are compact with Type I or II 
morphology, while all galaxies in the lower panel have type II morphology. Images are 3" to a side, and oriented with North up and East 
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set to 27.3 AB pixel"! (21.8 AB arcsec ). Red lines in the lower panel indicate the width and position angle of the LRIS-B slit during 
observations of each galaxy. 



10 



Law et al. 



TABLE 3 

Correlations of Bulk Outflow Velocity and 
Galaxy Properties'* 







Std. Dev. 






0.6 






1.6 






1.5 


Age'' 




1.9 


SFRf 




-0.9 


SsFR^ 




-2.1 


E(B - 




0.7 


G 




-1.7 


M20 




1.6 


C 




-1.3 


A 




1.0 






1.7 






2.6 


rj 




2.3 


n 




-1.6 


b/a^ 




-0.3 


H\m 




0.3 



^ All values arc the number of standard deviations from 
the null hypothesis that the quantities are uneorrelated, 
based on a Spearman rank correlation test. Positive (neg- 
ative) values indieate positive (negative) correlation be- 
tween the quantities (i.e., larger r corresponds to weaker 
outflows (less negative A^)Is). The sample size is 35 galax- 
ies. 

^ Cold gas mass, estimated from the star formation sur- 
face density and the observed Ha size, as in Erb et al. 
(2006c). 

Stellar mass, estimated from SED fitting. 
Total baryonic mass Mbar = M, -|- Mgas . 
Average population age from a constant SFR SED 
model. 

^ Star formation rate measured from SED fitting. 

^ Average star formation rate surface density Esfr. — 

SFB./Tvrl. 

^ Color excess due to dust extinction from a constant SFR 
SED model. 

^ GALFIT semi-major axis radius. 

J GALFIT circularized effective radius. 

^ Morphological minor/major axis ratio. 

detail we focus on the HAHQ sample for which accurate 
systemic redshifts have been obtained from rest-optical 
nebular emission. The 35 galaxies in our HAHQ sample 
represent t he subset of t he 73 galaxies previously con- 
sidered by ISteidel et alJ (pOinh for which 775T/WFC3 
rest-optical imaging data has now been obtained. In 
their larger sample, ISteidel et alJ (|201Clf ) demonstrated 
that there is a correlation between baryonic galaxy mass 
Mbary and Auig, the mean velocity offset of the inter- 
stellar absorption lines (zig) from the systemic redshift 

{zHa)- 

Awis = c(zis - Zsys)/ (1 + ^^sys) (1) 

Using th e Spe arman rank correlation test, 
ISteidel et al.l ()2010[ ) found that the null hypothesis 
(that Az;is and Mbary are uneorrelated) could be 
ruled out at 2.7 a signifi c ance. We repeat the cal- 
culations of ISteidel et alJ (|201Cll ) using an estimate of 
Mbary = M^-I-Mgas, where is the stellar mass derived 
in Paper I and Mgas the gas mass estimated using NIR- 
SPEC and/or OSIRIS H a fluxes in combin ation with the 
Ke nnicutt-Schmidt law ()Kennicuttlll998at see discussion 
bv lErb et aLll2006cD . In Tabled we list the confidence 
at which the null hypothesis for zero correlation can 



be rejected for Auig vs a variety of morphological and 
other physical properties. As expected w e recover the 
sense of the previous correlation noted by ISteidel et al.l 
()2010t ) between Auig and total baryonic mass, although 
this correlation has a significance of only 1.5cr; random 
resampling tests indicate that this lower significance is 
consistent with the lower total number of galaxies in our 
sample than in Steidel et al. (2010). 

4.1. Correlation of Low-ionization Absorption Line 
Strength with Semi-Major Axis Radius 

We find that Avis correlates with the semi-major axis 
radius r with a statistical significance of 2.6a in the sense 
that the largest velocity difii'erences (i.e., most negative 
At;is) are observed for the smallest (Type I) galaxies, 
while the largest (Type HI) galaxies have mean gas- 
phase velocities consistent with km . Combin- 
ing measurements from the individual low-ionization ab- 
sorption lines (except FeiiA1608, which is significantly 
weaker and noisier than the other transitions) we find 
that (Auis) = —226 ± 49 km s^^for the small radius 
half of the HAHQ sample (r < 1.95 kpc) compared to 
(Avis) = —125 ± 25 km s^^for the large radius half 
(r > 1.95 kpc)I3 We illustrate this trend visually in 
Figure 13 when sorted according to Awig galaxies with 
the smallest (blue boxes) and largest (red boxes) sizes 
are clearly distinguished. 

Adopting uncertainties of ^ 130 km s^^ in Auig and 
^0.1 kpc in r for each galaxy (Steidel et al. 2010; Pa- 
per I), we use a least squares routine to fit a first-order 
polynomial to the data and obtain the average relation 

Avis = (77-^- 342)kms"i (2) 
kpc 

We stress, however, that Avis does not characterize the 
outflow velocity of a particular galaxy, but rather repre- 
sents the centroid of a single-component fit to interstel- 
lar absorption lines that are a complex blend of mul- 
tiple components. In Figure [3] we show the stacked, 
continuum-normalized spectra for the small- and large- 
radius subsamples (solid blue/red lines respectively) in 
zoomed-in regions around the major low-ionization gas 
absorption lines. Careful inspection of Figure [9] sug- 
gests that the low ionization absorption lines in both 
low-r and high-r stacks have similar optical depths in 
their blue wings, indicating that both have outflowing 
gas tracing a similar range of velocities. Rather, the ma- 
jor difference is the growth of an absorption component 
centered at u km s~^ in the large-radius sample that 
shifts the centroid of the absorption line profile closer to 
the systemic velocity (see particularly AliiA1671) and 
increases the depth of composite absorption line pro- 
files even when precise systemic redshifts are unknown 
(e.g.. Figure lU top panel). This finding is similar to 
that described by Steidel et al. (2010; see also Weiner 
et al. 2009) who found that galaxies below a threshold 
Mbary « 4 X 10^°Mq lack the w ~ km s~^ component 
to the low-ionization interstellar absorption lines present 
in higher-mass galaxies. 

Following the formalism adopted by Steidel et al. 
(2010), we define the additional optical depth At{v) that 

^•^ Uncertainties represent la uncertainties in the mean combin- 
ing A?;is of individual galaxy spectra. 
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Fig. 7. — Mean rest-frame velocity offset AtJjs between the systemic (Ha) redshift ZHa and the low-ionization interstellar absorption line 
redshift zis as a function of major-axis galaxy size r (left panel), baryonic mass (middle panel), and SFR surface density EgpR (right panel) 
for the subset of 35 galaxies with systemic redshifts from Ho observations and high quality UV spectra. The arrow indicates the location 
of Q2343-BX453 {Avis = -944km s"! ) which along with Q2343-BX587 lies in the region for which r is unresolved at the 3cr level (grey 
shaded region). Crosses, open boxes, and open circles represent galaxies with morphological Type I, II, and III respectively. The vertical 
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would need to be added to the low-r composite galaxy 
spectrum in order to produce an absorption line profile 
identical to the high-r composite spectrum: 



/hr(l')=/lr(l')e-^^('') 



(3) 



where Ihriv) and Iir{v) are the spectral intensity of the 
high-r and low-r halves of the HAHQ sample respec- 
tively. In Figure[TU]we show At{v) for the three strongest 
unblended low-ionization absorption lines (SiiiA1260, 
C II A1334, and Si ii A1526) along with the average com- 
posite for all three lines combined. Although there is 
some variation between the individual lines, the com- 
posite exhibits a significant broad absorption component 
from ^ —400 to ~ -1-400 km s~^ in galaxies with the 
largest semi-major axis radii. As discussed at length by 
Steidel et al. (2010), this component could represent ei- 
ther infalling gas or stalled winds at large radii, although 
it most likely corresponds to cold gas located at small 
radii within the galaxy. 

4.2. Internal Correlation between the Morphological 
Parameters 

As indicated in Figure [TT] there are significant under- 
lying correlations between r and the other morphological 
statistics. On average, galaxies with larger semi-major 
axis radii have larger M,, , shallower radial density profiles 
(lower n), less concentrated light profiles (i.e., lower C 
and G), and tend to be more irregular (whether measured 
using M20, A, or ^>). Such generic trends are apparent 
in Figure [SJ the largest galaxies tend to be those that 
show most morphological irregularity to the light pro- 
files, while the smallest tend to be single well- nucleated 
sources. 



The physical mechanism underlying the change in Avis 
is the buildup of neutral gas at rest with respect to 
the stars (as discussed further in fjl])- The internal cor- 
relation of semi-major axis radius with other parame- 
ters however gives rise to a variety of correlations of 
marginal significance (see Table [3]) between Avig and 
stellar/baryonic mass, the six morphological statistics 
G, M20, C, A, 5", r, and n, and the SFR surface den- 
sity SsFR- Given the known positive correlation between 
mass and effective radius for the galaxies in our sample 
(Nagy et al. 2011; Paper I), it is perhaps unsurprising 
that similar qualitative trends in UV spectral features are 
observed when the galaxy sample is subdivided accord- 
ing to either mass or radius. Likewise, although there 
is no significant correlation of Avis with the total SED- 
derived star formation rate (in part since all galaxies in 
our sample are selected on the basis of rest-UV color to 
be actively star forming with SFR > 2Mq yr~^), there 
is a correlation at the > 2a level with the SFR surface 
density in the sense that galaxies with higher Ssfr tend 
to be missing the w ~ km s~^ absorption component 
(see Figure [71 right-hand panel) 13 Similarly, there is a 
^ 2(7 correlation of Auig with stellar population age in 
the sense that older populations have the strongest v ^ 
km s~^ absorption component), consistent with our find- 
ing that the r < 1.95 and r > 1.95 kpc samples have 
mean ages of 260 and 740 Myr respectively. 

4.3. High-ionization absorption lines 



Kornei et al. (2012) find that Sgp^ is more strongly correlated 
than galaxy size with Avis for ^ sample of star forming galaxies at 
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Fig. 8. — -ff5T/WFC3 F160W morphologies of galaxies in the HAHQ sample (i.e., with Ha-derived systemic redshifts and high-quality 
UV spectra) sorted (left-right, and top-bottom) in order of declining velocity offset Atijg. Blue boxes denote the three galaxies with the 
smallest semi-major axis radii in the HAFIQ sample (Q1623-BX453, Q2343-BX587, Q1217-BX95), red boxes denote the three galaxies with 
the largest semi-major axis radii in the sample (Q2343-MD59, Q2343-BX442, Q2343-BX389). The largest three galaxies all have formal 
AiJjg > km s~^ , while the smallest tend to be those with the most negative —950 < A-uig < —300 km s~^ . As demonstrated in Figurel4l 
the small-r sample has a 'younger' UV spectrum with significant Lya emission. Individual postage stamps are 3" to a side, oriented with 
North up and East loft. 



In addition to the low-ionization absorption lines, the 
high-ionization transitions SiivA1393, SiivA1402, and 
CivA1550 are also systematically less blueshifted for 
larger galaxies. Combining measurements made from 
the two Si IV absorption linesf^ we find that the mean 
blueshift of the high ionization lines (HIS) is (Awhis) = 
— 147 ± 58 km s~^ for the small-radius sample, but only 
(Awhis) = — 21±40 km s^^ for the large-radius sample 
In contrast to the low-ionization lines however, this effect 
is not produced by an additional component of t centered 

We neglect C IV in this stack since it is complicated by the 
superposition of interstellar absorption on a stellar wind-induced 
P-Cygni profile. 

Uncertainties represent Icr uncertainties in the mean combin- 
ing Avjjis of individual galaxy spectra. 



around w km s~^ , suggesting that the w '~ km s~^ 
gas may be more strongly self-shielded than the majority 
of the outflowing gas. Instead, the large-radius subsam- 
ple exhibits a weaker blue wing to the absorption proflle 
(Figure ITUl) . suggesting that highly ionized (T > 10* K) 
gas may have lower peak outflow velocities in galaxies 
with rest-optical continuum radii r > 3 kpc. 

4.4. Fine structure emission 

In addition to exhibiting stronger Lya emission and 
more systematically blueshifted low-ionization absorp- 
tion features, galaxies with smaller rest-optical contin- 
uum radii also have stronger flne structure emission lines. 
As indicated by Figure |1 Si ii* A1265 and A1309 emis- 
sion features are more prominent in composite spectra of 
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Fig. 9. — As Figure |4] but zoomed in to the central 1500 km s ^ surrounding each of the major low-ionization 
(Si II A1260, O I A1302+Si ii A1304, C ii A1334, Si ii A1526, Fe ii A1608, Al ii A1670) and high-ionization (Si iv A1393, Si iv A1402, C iv A1550) ab- 
sorption Une features. Blue/red Hues respectively represent the stacked spectra of galaxies in the smallest/largest size bins. Shaded regions 
indicate la uncertainties estimated by Monte Carlo bootstrap resampling from the individual galaxies making up the composite spectra. 



galaxies with effective radii r <2 kpc than in those with 
r > 3 kpc (most noticeably in the composite stack from 
the parent sample, and to a lesser extent in the HAHQ 
composite stack). This result is consistent with the re- 
cent study of Jones et al. (2012), who noted that Sill* is 
stronger in star forming galaxies at z 4 than at lower 
redshifts and ascribed decreasing Sill* equivalent width 
to an increasing characteristic radius of low ionization 
gas. 

5. DISCUSSION 

ISteidel et al.l (|2010( ) suggested that the physical prop- 
erty governing both Lya emission properties and inter- 
stellar absorption line kinematics was baryonic mass, 
which correlated with A wis with a significance of 2.7 a 
in their sample of 73 galaxies. Since galaxies with higher 
baryonic masses tend to have older stellar populations 
and greater total quantities of interstellar gas, this ad- 



ditional gas may naturally lead to greater covering frac- 
tions and longer scattering paths, suppressing Lya emis- 
sion and making interstellar absorption lines deeper. 

In detail though, mass may not be the most predic- 
tive indicator of the spectroscopic properties of a given 
galaxy. Although the correlation that we find between 
Avis E^nd r is of comparable statistical significance (2.6cr) 
to the baryonic mass correlation described by Steidel et 
al. (2010), this may simply be because there are ~ a 
factor of 2 fewer galaxies in the morphological HAHQ 
sample. Replicating the HAHQ sample to a total sam- 
ple size of 73 galaxies with random resampling using 
a Monte Carlo algorithm (repeated 10,000 times), we 
estimate that the expected correlation significance be- 
tween r and Ai;is in a sample of the size of Steidel 
et al. (2010) is 3.8 ± 0.7 ct. Indeed, while the galax- 
ies with the strongest/ weakest Avis occupy an overlap- 
ping range of Afbary (the galaxy with the weakest Avis 
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Fig. 10. — Additional optical depth At{v) required to transform the low-r composite galaxy spectrum into the high-r composite as 
a function of offset from the systemic velocity. Left-hand panel: Colored lines indicate the profiles derived for individual low-ionization 
interstellar (LIS) absorption lines, solid black line indicates the average of the three individual At{v) profiles. Right-hand panel: Colored 
lines indicate the profiles derived for SilvA1393 and A1402 high-ionization interstellar (HIS) absorption lines, solid black line indicates the 
average of the two individual At{v) profiles. Note that the LIS profiles show additional optical depth near v ~ - 200 km s~l in large 
galaxies, while the HIS profiles show marginally reduced optical depth at relative velocities —600 < v < km s~^ . Shaded grey regions 
indicate la uncertainties estimated by Monte Carlo bootstrap resampling from the individual galaxies making up the composite spectra. 



in our HAHQ sample has a baryonic mass smaller than 
the galaxy with the strongest Avis), these galaxies have 
clearly distinct visual morphologies (Figure [8]). Like- 
wise, although the galaxies with the highest Ssfr are 
those with the strongest Ai;ig (consistent with theoret- 
ical expectations that galaxies with a higher concentra- 
tion of supernovae should be capable of more completely 
expelling the local ISM), this correlation is driven en- 
tirely by the (stronger) relation between r and Awjg and 
the definition Ssfr = SFR/r^ (although c.f. Kornei et 
al. 2012 for star forming galaxies at z ~ 1). 

All of the correlations detailed herein may therefore 
simply be byproducts of a single underlying physical 
mechanism; the growth of an extended disk. Star- 
forming z ^ 2 — 3 galaxies with low to moderate stellar 
mass (i.e., < lO^^M©) tend to be compact, young, 
triaxial systems whose kinematics are dominated by ve- 
locity dispersion rather than rotation about a preferred 
kinematic axis. High SgpR within these systems drives 
strong gaseous outflows into the surrounding IGM and 
suppresses the buildup of a stable ISM. In time, the in- 
creasing stellar mass formed by these galaxies is able to 
stabilize the formation of an extended (albeit still thick) 
gaseous disk (as proposed by, e.g., Martig & Bournaud 
2010) in which the overall gas fraction decreases (/gas 
is anticorrelated with Awig with a confidence of 2.8a in 
the HAHQ sample) and systemic rotation becomes an in- 
creasingly important means of dynamical support. The 
lower SsFR in this extended disk is no longer able to ex- 
pel gas as efficiently, a large component of which remains 
close to the galaxy and serves to both extinguish the Lya 
emission (or resonantly scatter it to larger radii) and su- 
perimpose a zero- velocity absorption component atop the 
stellar continuum emission. 

While the growth of such extended structure is natu- 
rally associated with both evolving and Ssfr, there 



is not a single well-defined value of either that demarcates 
the transition, and the mass at which extended structures 
may be supported can vary significantly based on the 
merger histories and other details of individual galaxies. 
In contrast, rest-frame optical size r characterizes where 
the transition to an extended disk has already occurred. 
From a phenomenological standpoint, rest-optical mor- 
phology may therefore be the most robust and accessible 
means of identifying systems for which feedback is less 
efficient at evacuating gas from the galaxy. 

This physical picture necessarily suggests that the 
zero-velocity gas feature should become more pro- 
nounced (i.e., Avis is less strongly blueshifted) as the 
contribution of rotational shear to the dynamical sup- 
port of a galaxy increases. Such trends have already been 
noted in samples of galaxies observed with either long-slit 
(see Figure 2 of Steidel et al. 2010) or IFU spectroscopy 
(see Figure 13 of Law et al. 2009). More recently, Forster 
Schreiber et al. (2011a) presented HST/mCMOS mor- 
phologies of a sample of six galaxies observed previously 
with the SINFONI integral field spectrograph. Of these 
six galaxies, five are representative of the most disk-like 
known systems at z ~ 2 based on their observed kine- 
matics (the other galaxy is identified as a major merger). 
These five galaxie^ all lie at the upper edge of the 
distribution of effective radii for z ~ 2 — 3 star form- 
ing galaxies (see colored circles in Figure 13 of Paper I) 
and have a mean (A^;is) = -4 km s'^ (excluding Q1623- 
BX663, whose rest-UV spectra are complicated by AGN 
features). These systems can therefore be identified from 
the overall z ~ 2 — 3 galaxy sample by both their semi- 
major axis radii and spatially resolved velocity shear, but 

I'' Q1623-BX663, SSA22a-MD41, Q2343-BX389, Q2343-BX610, 
and Q2346-BX482. One of these galaxies, Q2343-BX389, is also in 
our WFC3 sample. 



Relation ol z ^ 2 Morphologies to Gas Kinematics 



15 



4 - 



o 




2 4 
C 



0.2 0.4 0.6 0.8 2 4 
A ^ 



12 3 
2= 



■■SFR 




O 





9 9.5 10 10.5 0.3 0.4 
log(M,/Mj G 



20 



1 2 3 
n 



Fig. 11. — Correlations between semi-major axis radius r, stellar mass M«, and the other morphological parameters. Symbols are as in 
Figure [T] Each panel quotes the statistical significance of the correlation on the basis of the Spearman rank correlation test. 



not by their stellar masses, which are widely distributed 
in the range = 8 x 10^ - 1 x lO^^M©. 

Steidel et al. (2010) suggested that the anticorrela- 
tion between Avig and velocity shear may be due in part 
to projection elFects. If objects with the largest veloc- 
ity gradients tend to be edge-on rotating disks, outflows 
from these galaxies might be expected to be coUimated 
perpendicular to the disk so that there is simply little 
projection of the outflowing gas onto our line of sight 
(see, e.g., Chen et al. 2010, Bordoloi et al. 2011), while 
face-on systems would exhibit minimal rotation and max- 
imum Awis- Such a correlation between inclination and 
outflow velocity was recently found for a sample of star 
forming galaxies at z ~ 1 by Kornei et al. (2012). How- 
ever, in our z ~ 2 sample we observe no correlation of 
Avis with inclination as parametrized through the appar- 
ent minor/major axis ratio b/a for the 35 galaxies in our 
HAHQ sample (0.3cr from the null hypothesis; see Table 
[3]). In part, the lack of correlation may simple reflect the 
irregularity of z 2 star forming galaxies and the difS- 



culty of determining inclination robustly for these intrin- 
sically triaxial (Paper I) systems. However, no significant 
difference is observed in Avis even for the two most disk- 
like galaxies in the sample. Morphology and IFU-derived 
kinematics suggest that Q2343-BX389 is consistent with 
an edge-on disk (Forster Schreiber et al. 2009), while 
Q2343-BX442 has both morphology and IFU kinemat- 
ics consistent with a nearly face-on (i = 42° ± 10°) disk 
with spiral substructure (see discussion by Law et al. 
2012b). Despite their different orientations both galax- 
ies have Avis ~ km , suggesting either that Avis is 
uncorrelated with inclination, or that the opening angle 
of the outflow is less than ~ 40°. 

As illustrated by Figure [10] much of the additional 
absorption in the largest galaxies has positive velocity 
'--^ — 400 km s~^ with respect to the systemic redshift, 
and formally Avis > km s^^ for the three galaxies 
with largest semi-major axis radii in our HAHQ sam- 
ple (Q2343-MD59, Q2343-BX442, Q2343-BX389). It is 
therefore possible that absorption line spectroscopy of 
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the largest galaxies is tracing infalling gas, whether from 
a recent merger, recycled wind accretion (e.g., Oppen- 
heimer et al. 2010), or from cold cosmological flows (e.g., 
Dekel et al. 2009; Keres et al. 2009). While such inflows 
are one plausible mechanism for solving the gas fueling 
probem, it is unclear why such accretion (presumably 
with a small covering fraction) would cause the near- 
complete absorption of Lya (see discussion by Steidel et 
al. 2010). Similarly, simulations (e.g., Dekel & Birnboim 
2006; Keres et al. 2009) generally predict that such cold 
flows should be suppressed in the most massive galaxies 
at z > 2, in contrast to our observation that Avis > 
km occurs more commonly for large and massive 
galaxies. 

6. SUMMARY 

We have demonstrated that rest-optical morphology 
is correlated with the gas-phase properties of z ~ 2 
star forming galaxies as traced by their rest-UV spec- 
tra. Lya emission is most commonly observed in galax- 
ies with small rest-optical half-light radii (r < 2 kpc), 
and the equivalent width W^Lya of the emission is corre- 
lated with half light radius with > 2<t confidence. Al- 
though some large-r galaxies with Type III morphology 
(i.e., extended and diffuse) show Lya in emission, this 
emission is relatively weak; the twelve galaxies in our 
sample with VFlyq > 30 A all have small half-light radii 
(r < 2 kpc) and Type I or Type II morphologies (i.e., 
compact clumps) . This correlation is independent of the 
better- known relations between VFlyq, dust content, and 
star formation rate. Additionally, all galaxies in our spec- 
troscopic sample for which Lya is observed in emission 
and no absorption lines are present in the rest-UV spec- 
trum have morphological Type II (i.e., consist of two or 
more clumps of emission separated from each other by ^ 
1 arcsec), suggesting that slit losses of continuum light 
from these galaxies may be significant. 

Combining our res ults with those previously published 
in the liter ature Ce.g..[Steidel et al.ll2010l iPentericci et al.l 
[2010, Kornd_eLalJl2Q10D, the broad picture of z ~ 2 - 3 
star forming galaxies is relatively clear. Galaxies with 
Lya emission tend to be physically smaller than galaxies 
with negligible Lya emission, with lower stellar mass, less 
dust, marginally younger stellar populations, lower star 
formation rates, and a slightly greater degree of nucle- 
ation in the rest-frame UV and optical light distribution. 

Using a subsample of 35 galaxies with both high- 
quality rest-UV spectra and Ha-derived systemic red- 
shifts we have found that galaxy morphology is also 
correlated with the median blueshift of interstellar ab- 
sorption line features Awis with 2.6(7 significance. Al- 
though similar blue wings in the low-ionization absorp- 
tion line profiles indicate that galaxies of all morpholog- 
ical types drive comparably strong outfiows at speeds 
up to ~ 800 km s~^ into their surrounding environ- 
ments, increasing optical half-light radius is accompa- 
nied by an increase in optical depth around the systemic 
redshift of the galaxy. This finding is similar in statisti- 
cal significance to tre nds with incre asing baryonic mass 
described recently bv iSteidel et al.l (2010), but may be 
of physically greater importance considering the smaller 
sample size. While galaxies with Avis differing by al- 
most 1000 km s""'^ can have baryonic masses identical 
to within observational uncertainty, such galaxies have 



clearly distinct morphologies spanning the range from 
compact clumps to extended disks. 

The increased optical depth arises in gas nearly at rest 
with respect to the systemic redshift, in close physical 
proximity to the galaxy, and apparently corresponding 
to an extended disk of optical continuum light. This 
stronger w ~ km component is responsible both 
for the shifting absorption line profile and the attenu- 
ation of Wi,ya (via a longer resonant scattering path) 
in galaxies with larger rest-optical half light radii. The 
V ^ km s~^ absorption component is not present in 
the higher-ionization features SiivA1393 and 1402 how- 
ever, suggesting that this gas may be more strongly self- 
shielded than the majority of the outflowing gas. 

In contrast to r ecent results at lower redshifts (e.g., 
IKornei et al.ll2012l ). there is no obvious relation between 
Awis and galaxy inclination as parameterized by the axis 
ratio b/a for the 35 galaxies in our HAHQ sample. In 
part this may be due to the difficulty determining ro- 
bust inclinations for inherently triaxial star-forming sys- 
tems. Even for two particularly well-resolved galaxies 
with visual morphology and IFU kinematics consistent 
with face-on and edge-on disks however, there are no 
obvious differences in the outflowing gas kinematics (al- 
though the quality of individual spectra prohibits a de- 
tailed analysis of just these two galaxies). There is there- 
fore no evidence to support a classical bipolar outflow 
model for star forming galaxies at z ~ 2 — 3. This may be 
because such galaxies are in the process of transforming 
from clumpy, irregular systems to more regular (albeit 
thick) disks, and outflowing gas at high velocities and 
large distances from the galaxy may predate establish- 
ment of the disk and its corresponding polar axis. 

In general, our observations are consistent with inside- 
out growth of star forming galaxies in the young uni- 
verse: 'Typical' z ^ 2 star forming galaxies appear to 
be gas-rich, compact, triaxial systems systems that are 
dominated by velocity dispersion between individual star 
forming regions rather than systemic rotation, and whose 
high SsFR drives strong outflows into the surrounding 
IGM. As these galaxies mature they gain stellar mass, 
stabilizing the formation of extended (albeit still thick) 
gaseous disks in which rotational support is observed to 
play an increasing role. As the star formation migrates 
from central regions into these extended disks the Ssfr 
drops and is no longer capable of expelling gas from the 
galaxy as efficiently. 

We caution, however, that it is challenging to gener- 
alize about the relation between rest-optical morphology 
and gas phase kinematics for the entire z ~ 2 star form- 
ing sample on the basis of the relatively small number of 
galaxies in our sample. A major limiting factor is sys- 
temic nebular redshifts, without which we cannot quan- 
tify the outflow velocities. Although ~ 100 z ~ 2 — 3 
star forming galaxies with rest-UV spectra in the KBSS 
have had systemic redshifts derived from nebular emis- 
sion line spectroscopy (e.g., Erb et al. 2006b; Law et 
al. 2009), only a fraction of these galaxies fell within the 
footprint of our HST/WFC3 imaging survey. In the fu- 
ture it will be possible to explore these trends in much 
greater detail with the aid of multi-object near-IR and 
optical spectrographs, using which it will be possible to 
obtain Ha redshifts and outflow kinematics efficiently for 
large and diverse samples of galaxies in fields that have 
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already been imaged with HST. 
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